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Abstract. We report on multi-epoch HSTAVFPC2 images of the XZ Tauri binary, and its outflow, covering the period from 
1995 to 2001. Data from 1995 to 1998 have already been published in the literature. Additional images, from 1999, 2000 and 
2001 are presented here. These reveal not only further dynamical and morphological evolution of the XZ Tauri outflow but 
also that the suspected outflow source, XZ Tauri North has flared in EXor-type fashion. In particular our proper motion studies 
suggests that the recently discovered bubble-like shock, driven by the the XZ Tauri outflow, is slowing down (its tangential 
velocity decreasing from 146 km s^^ to 1 17 km s^^). We also present simulations of the outflow itself, with plausible ambient 
and outflow parameters, that appear to reproduce not only the dynamical evolution of the flow, but also its shape and emission 
line luminosity. 
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1. Introduction 

XZ Tau is a classical T- Tauri binary system with a separation of 
(X'S JHaas et al 19901 1 and is located in the well known Lynds 
1551 star-forming region some 140pc away (Elias 1978 1. The 
system was first found to have an associated Herbig-Haro 
(HH) outflow through ground-based CCD imaging and spec- 
troscopy ( Mundt et al 1988 , Mundt et al 1990). These early ob- 
servations revealed a bipolar optical flow that could be traced 
to at least 10" on either side of the binary at a position angle 
of 15°. The first Hubble Space Telescope (HST) Wide Field 
Planetary Camera 2 (WFPC2) images of XZ Tauri taken in 
1995 show a bubble of emission nebulosity extending 4" to the 
north of the system (Kristet al 1997 hereafter K97). Further 
images, taken 3 years later, show dramatic structural changes as 
the bubble expanded and altered from being centre-brightened 
to limb-brightened, suggesting the formation of a HH bow- 
shock jKrist et al 1^9^ hereafter K99). 

Ground-based photometry of XZ Tau from 1962 to 1981 
JHerbst et al 1994 ) has shown variations, of almost two magni- 
tudes in the V band, for the binary as a whole. Such variations 
are common amongst young stellar objects (YSOs). Of the two 
components, the southern one has been, at least until recently, 
optically brighter and thought to be the more evolved star. Its 
companion, however, dominates at infra-red wavelengths and is 
probably of higher overall luminosity (Haas et al 1990). Recent 
Faint Object Spectrograph (FOS) observations, however, unex- 
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pectedly found the northern component to be optically brighter 
J White & Ghez2001 hereafter WGOl), a result that we will 
discuss further in the light of our findings. For this system 
therefore it seems more appropriate not to use the terms pri- 
mary and secondary but instead we will adopt the nomencla- 
ture, used in K97, of XZ Tau North and South. 

We report here analysis of further HST Archive WFPC2 
images of the XZ Tau system and outflow from 1999, 2000 and 
2001 . These data show not only ongoing changes in the outflow 
but a dramatic brightening of XZ Tau North in the optical sug- 
gesting that it may be an EXor. We also simulate the outflow in 
an attempt to reproduce its dynamical and morphological evo- 
lution. 



2. Data 

High resolution archival WFPC2 images of XZ Tauri were 
obtained for 3 consecutive epochs: 1999 February 3; 2000 
February 6; and 2001 February 10. The binary was at the 
same approximate location on the Planetary Camera (spatial 
sampling = 0'.' 04555 pixel" ^) for all frames. Four filters were 
used: Ha (F656N); [SII] (F673N); R-band (F675W) and I-band 
(F814W). The frequency and duration of the short exposures 
for each filter were: lxl20s; lxl80s; 2x6s; and 2x6s re- 
spectively. Two long exposures for each filter of 1000 s were 
also extracted from the HST Archive, all of which were sat- 
urated at the location of XZ Tauri. No short exposures were 
made for the I-band filter in the final year. 
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Previously published WFPC2 archival data (K99) for 
XZ Tauri facilitated proper motion measurements. These data 
comprised images from 1995 January 5 and 1998 March 25. 
The 1995 images used were 2 x 600 s exposures in the R-band 
and 1 X 3 s exposure in the I-band. The latter allowed us to de- 
termine the stellar positions, as no R-band short exposures were 
taken and the long exposures were saturated in the vicinity of 
the star. The 1998 images were taken in the R-band only: 2 x6 s 
and 2 X 1100 s exposures. All frames were processed through 
the standard HST pipeline and each set of double exposures 
was combined to eliminate cosmic rays. 

The accumulated data yielded high resolution images of 
XZ Tau spanning a total of 6 years. We analyse these data 
in Section 3, under the two headings: Outflow Structure, 
Proper Motions & [SII] Luminosity and Stellar Astrometry & 
Photometry. 

3. Results 

3. 1. Outflow Structure, Proper Motions k [SII] 
Luminosity 

The Planetary Camera images for each year were aligned and 
their orientation corrected where necessary. A montage of R- 
band images, covering the full timespan, is shown in Fig. ^ 
As previously observed (K99), there appears to be two edge- 
brightened "bubbles", i.e. inner and outer shocks, similar to 
those seen with HST close to DG Tau ( Bacciotti et al 2000) . 
By 2001, the rapidly fading outer shock reached a distance of 
approximately 800 AU from the binary, while the inner shock 
has travelled approximately half as far 

Although remarkably bright in 1995, the knots in the 
XZ Tau jet appear very faint by 2001. The jet's PA is around 
15 ° in line with the major axis of the elongated outer bubble 
and the known extended optical outflow. The knots however are 
not in a perfectly straight line as some deviation is evident and 
reminiscent of similar behaviour seen in flows like the HH 34 
jet close to its source p^^eipurth et al 200*2) . 

Significant proper motions were observed as the bubble ex- 
panded away from the XZ Tau. Tangential velocities for the 
central jet knot and outer shock (see for example the bottom 
right "difference" frame in Figure Q were obtained from the 
long F675W exposures. The latter were used since they in- 
clude the main HH emission Unes, i.e. [OI]AA6300,6363, Ha 
and [SII] A A67 17, 6731. The distance to the binary was assumed 
to be 140 pc as in K99. In Table^we list the derived tangential 
velocities as a function of epoch for the outer shock and an av- 
erage velocity for the jet. Longitudinal sizes and speeds were 
calculated along the jet axis (P.A.~15 ° ). Transverse sizes and 
speeds were calculated perpendicular to the jet axis at the bub- 
ble's widest point, (specifically, the point where the source-to- 
shock line forms an angle of 20 ° with the jet axis in the 1995 
image). As the shock expands it grows fainter and so its width 
is not given for later observations. The central knot also grows 
very faint, especially at later epochs, so we prefer to quote an 
average for the jet speed. Finally, although the inner shock or 
bubble was seen to expand, changes from year to year proved 
difficult to trace and so no tangential velocity is given here. 



Observation 


Distance 


Bubble 


Longitudinal 


Transverse 


Date/Interval 


from source 


Width 


Speed 


Speed 






f AT 11 


(kms~^) 


(km s~ ^) 


1995 Jan 5 


598 ±16 


336 ±25 






1995 Jan 5 - 1998 Mar 25 


697 ±10 


405 ±14 


146 ±28 


101 ±42 


1998 Mar 25 - 1999 Feb 3 


721 ±9 


423 ±5 


130 ±73 


105 ±83 


1999 Feb 3 -2000 Feb 6 


746 ±11 




119 ±65 




2000 Feb 6-2001 Feb 10 


771 ±10 




117 ±69 




1998 Mar 25 - 2001 Feb 10 






130 ±84 (jet) 





Table 1. Projected sizes and speeds of the XZ Tauri outer bow- 
shock (and an average speed of the jet knot) as it evolves. 
Distances are quoted for the latter date of the observation in- 
terval. Errors are at the three sigma level. 



From Table ^ it is clear that we see a deceleration of 
the outer shock. Although the errors appear large, if we con- 
sider the interval 1998-2001 we obtain a longitudinal speed 
of 121 (± 24)kms-i (cf. 146(± 28)kms-i from 1995 to 
1998) suggesting the shock front is decelerating. Simulations 
were carried out in order to model the observations using plau- 
sible jet and ambient medium parameters (see §3.3 below). 
Observerd speeds are as expected in order to maintain the ob- 
served low aspect ratio of the bubble, and are in the same re- 
gion as those given by the simulation, (see Table0]below). As 
a further check the simulated and observed [SII]AA6317,6731 
luminosities were compared. The latter was found from the 
2001 F631N images to be 2.5 (± 3.6) x lO^s erg s^^, exclud- 
ing the inner shock region which is contaminated with diffrac- 
tion spikes and nebulosity near the star If we assume an ex- 
tinction of A„~1.39 towards the bubble, i.e. the same as to- 
wards the binary (WGOl), the intrinsic luminosity increases to 
--4.3 xlO^^ergs-^ 

3.2. Stellar Astrometry k Photometry 

Using the short R band exposures, the separation and position 
angle of the binary were determined for the different epochs 
and are listed in Table |2] Within errors, no changes were de- 
tected in the separation of the XZ Tau binary although an av- 
erage decrease of 0.75 ° yr^^ was observed in its position an- 
gle, a value somewhat higher than derived by K99 based on 
earlier HST data alone (of 0.5 °yr^^) but at the same time 
lower than that found by |Woitas, Kohler, & Leinert (2001)| us- 
ing speckle interferometry (of ~1.3°yr~^). Following K99, 
if we assume a face-on orbit, the total mass for the binary 
is about 0.3 Mq. Although such a value is clearly too low, 
note that the combined mass is very sensitive to the projection 
angle. Such a figure is certainly lower than the estimates of 
Hai'ti gan'& Kenyon (2003)| and WGOl suggesting a combined 
mass closer to 1 M0. 

Broadband R and I magnitudes were calculated for the bi- 
nary, using the method outlined in the WFPC2 Data Reduction 
Handbook and interpolated Johnson offsets appropriate for the 
spectral types of XZ Tau North and South, i.e. M2 and M3.5 
respectively (Hartigan & Kenyon 2003). Our results are pre- 
sented in Table 13.21 and the R-band WFPC2 data were used 
to plot a light curve for each component. Fig. |2l Over the six 
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Fig. 1. HSTAVFPC2 F675W (R-band) images of the XZ Tau outflow on: 1995 Jan 5; 1998 Mar 25; 1999 Feb 3; 2000 Feb 6; and 
2001 Feb 10. The short exposure image of the binary is superimposed on the long exposure image in each case. The bottom right 
frame is an R-band difference image 1999-1998 showing proper motion in the jet and the outer shock. 



Date of 


Separation 


Position Angle 


Observation 


(arcsec) 


(degrees) 


1995 Jan 5 


0.302 


146.35 ±0.38 


1998 Mar 25 


0.299 


145.50 ±0.25 


1999 Feb 3 


0.303 


144.85 ±0.26 


2000 Feb 6 


0.300 


143.30 ±0.25 


2001 Feb 10 


0.298 


142.38 ±0.24 



Table 2. Separation and position angle of the XZ Tau binary. 
Errors in separation angles are estimated to be ±0.005. 



years of observations, XZ Tau South, although it varies, does 
so by at most 0.3 magnitudes (in R) and has a mean R magni- 
tude of around 13.5. In contrast, XZ Tau North shows an initial 
reduction in brightness of about a magnitude in R until 1998 
and thereafter its brightness increases dramatically by around 3 
magnitudes. This flaring behaviour means that by 2001 XZ Tau 
North was actually the brighter star Similar variations are seen 
in the I band although the data is somewhat more sparse. 

The dramatic brightening of XZ Tau North suggests it is 
an EXor. EXors, named after their prototype EX Lupi, are 
extreme classical T Tauri stars that periodically undergo out- 
bursts from the UV to the optical. Although increases by sev- 
eral magnitudes with rise times of up to a few years have 



been recorded [Herbig (1989)| the changes in these YSOs are 
not as extreme as in FU Orionis stars. EXor spectra, for ex- 
ample, even in outburst continue to resemble T Tauri stars. 
The phenomenon is thought to be due to major increases in 
the underlying disk accretion rate, but the number of EXors 
known is relatively small. The proposition that XZ Tau North is 
such a YSO is further strengthened by the HST spectroscopic 
data of |Hartigan & Kenyon (2003)1 As with other EXors the 
spectrum shows not only very strong Balmer lines but also 
strong Ca II, and moderate Na I in emission (see, for ex- 
ample, Parsamian, Mujica, & Corral 2002 1. Finally EXor be- 
haviour would explain why the the 1996 HST/FOS acqusi- 
tion of XZ Tauri (WGOl) unexpectedly locked onto the wrong 
(northern) star. WFPC2 images taken in March 1997, as part of 
HST GO Programme 6725, show XZ Tauri North comparable 
in brightness with XZ Tauri South at short optical wavelengths. 
This also suggests that EXor behaviour may be visible at the 
shortest wavelengths first. 



3.3. Modelling the Outflow 

In this section we present the results of a numerical simula- 
tion of the XZ Tauri outflow in order to check our physical 
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Filter 


Observation 


Flux Ratio 


XZ South 


XZ North 




Date 


Soiitli/North 


Magnitude 


Magnitude 


R-band 


1995 Jan 5 


5.11 


13.16 


14.93 




1997 Mar 8 


5.39 


13.47 


15.38 




1998 Mar 25 


9.02 


13.22 


15.68 




1999 Feb 3 


3.66 


13.45 


14.93 




2000 Feb 6 


2.43 


13.32 


14.36 




zOul Feb 10 


0.63 


13.25 


12.82 


I-band 


1995 Jan 5 


4.80 


12.03 


13.78 




1997 Mar 8 


6.30 


11.97 


14.02 




1999 Feb 3 


3.38 


12.04 


13.42 




2000 Feb 6 


1.98 


12.02 


12.82 



Table 3. Johnson apparent magnitudes of the XZ Tauri binary. 
EiTors in magnitude of ±0.05 were estimated based on the af- 
fect of changing aperture size, given that the stellar PSFs are 
overlapping. The 1995 Jan 5 R-band data is from K97 using 
Tiny Tim PSF fitting to the saturated R-band images. The 1997 
Mar 8 data is from GO Programme 6725. 
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Fig. 2. Stellar Johnson R-band apparent magnitudes of the 
XZ Tauri binary. 



interpretation of the observations. The code we use is that de- 
scribed in Downes & Ray (1999). It is a properly upwinded 
second order (in time and space), cylindrically symmetric code 
for solving the inviscid Euler equations. In addition to tracking 
the hydrodynamic variables (density, velocity and pressure), 
the code also tracks the ionisation state of hydrogen, without 
the assumption of ionisation equilibrium. It also has the capac- 
ity to track the number density of H2, but this number density 
was set to zero for these simulations. 

3.3.1 . Initial conditions 

The gas is taken to be one of solar abundances. The initial am- 
bient density is assumed to be 100 cm^^ close to the source, ris- 
ing to a value of 600 cm^'^ at a distance of 1.25 x 10^^ cm from 
the source, (see below for a discussion of why this behaviour is 
assumed). The ambient pressure was taken to be uniform, giv- 
ing a temperature of 10"^ K close to the source, and dropping to 
approximately 160K at the rise in ambient density. The resolu- 
tion was set to 1 X 10^^ cm in both the radial and poloidal (out- 
flow) directions (r and z respectively), and the total grid size 



was 450 X 1500. To estimate of the mass flux in the outflow, 
we assumed approximately 2% of the accreted mass ultimately 
ends up in the outflow ( Hartigan, E dwards, & Ghandour 1995) . 
A mass accretion rate for XZ Tau North of 1 x 10~^ M© yr~^ 
was used (Hartigan & Kenyon 2003) based on recent photo- 
spheric veiling measurements. 

In order to choose parameters for the outflow we must take 
the following issues into consideration: 

1 . The aspect ratio of the outer bubble/bowshock is small de- 
spite the fact that the outflow is thought to be highly su- 
personic with respect to the sound speed in the ambient 
medium, 

2. If we assume a constant jet/outflow radius based on the ob- 
served inner knots, i.e. r^ct ~ 2 x 10^^ cm, then the inferred 
mass-flux in the outflow would imply a large jet density and 
specifically one that is much higher than the ambient den- 
sity. This cannot be the case for the XZ Tauri outflow, as 
such a highly dense jet would simply plough through the 
local medium creating a very narrow bowshock, in contra- 
diction with our observations. 

One way of explaining the observations is to invoke a mod- 
erately coUimated wind. A wind with a moderate opening angle 
will have a density which decreases significantly with distance 
away from the source, and will also have a larger radius with 
distance. These effects each lead to a broader bowshock, in line 
with the observed morphology of the system. In addition, if the 
flow is not in the plane of the sky, the bowshock will appear to 
have a lower aspect ratio due to projection effects. We assume 
an angle of 30 ° between the plane of the sky and the axis of 
the outflow. 

With these considerations in mind, we set the full opening 
angle of the wind to be 22 ° with an initial diameter (full-width 
half-maximum in density) of 4 x lO^'' cm, or 40 grid cells. The 
wind density was given a linear profile across the outflow axis, 
with a density range from 1800cm~'^ along the jet axis to ^ 
100 cm^'^ at the edge, and was chosen so that the total mass- 
flux was 2 X lO~^Af0yr~^. Hence we are assuming the jet 
(traced by the knots) does not dominate the dynamics of the 
system. 

The observed tangential speed of the flow is in the range 
1 30-200 km s^^ so we choose a constant wind velocity of 
250 km s^^, which will give a tangential velocity for the jet 
fluid of roughly 216 km s^^. Since the bowshock appears to be 
decelerating significantly (see tabled, we need some way to 
slow it down. Turning off the jet will not, on its own, give us 
the observed deceleration since the momentum density of the 
material in the cooled shell between the bowshock and jetshock 
is very large. It was found that the most effective way of decel- 
erating the bowshock was to impose a significant increase in 
the ambient density. As mentioned above, a six-fold increase 
was found to roughly reproduce the observed velocities. 

3.3.2. Simulation results 

Figure |3l contains a greyscale plot of the distribution of den- 
sity at t — 21yrs, the time at which the simulation roughly 
matches the 2001 observations. The non-zero opening angle of 
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Fig. 3. Log-scale plot of the distribution of density at t ~ 
21yrs. The distance scales are in units of lO^'^ cm and the in- 
tensity scales are in units of g cm^^ 
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I 



note however that the bowshock apex appears much brighter 
with respect to the wings than is actually observed. This is most 
likely due to our use of cylindrical symmetry which results in 
a focussing of shocked material onto the axis in a way which 
would be very unlikely to happen in three dimensions. 

4. Conclusions 

Multi-epoch HST/WFPC observations of the XZ Tau binary 
and its associated outflow have shown considerable changes in 
the system within only 6 years, from 1995 to 2001. The pres- 
ence of two limb-brightened shock fronts is now clearly evi- 
dent, with a deceleration of the outer shock from 146 km s^^ 
to 117kms^^. Stellar photometry revealed that the suspected 
source of the outflow, XZ Tau North, has flared in EXor-type 
fashion increasing in brightness by 3 magnitudes in R be- 
tween 1998 and 2001. Finally, numerical simulations of the 
outflow produced reasonable agreement with observation in 
terms of morphology, dynamical evolution and emission line 
luminosity, using plausible ambient and outflow parameters. 
Deceleration by the amount observed, caused by the ambient 
medium, should have produced a much brighter bowshock apex 
than that seen. The cause of this discrepancy is not obvious. 



Fig. 4. Simulated [SII] image projected onto the sky (assum- 
ing an inclination to the plane of the sky of angle of 30°) 
at t = 21yrs (see text). The intensity scales are in units of 

erg s^^ cm^^ 



the jet is clearly noticeable. The bowshock itself has a number 
of irregularities which probably arise from the Vishniac insta- 
bility. Another possible source is the Rayleigh-Taylor instabil- 
ity, since the bowshock is decelerating. It can be seen that the 
bowshock has begun to encounter the rise in the ambient den- 
sity (see section 13. 3. 1> . This leads to a marked decrease in its 
speed of advance (see table |^. The behaviour of the velocity 
of the outer shock matches that observed reasonably well (see 
tablegj. 

Figure 0]is a simulated [SII] image calculated for the den- 
sity distribution shown in figure |3 assuming an angle to the 
plane of the sky of 30 °. In calculating this we assume that the 
ionisation fractions state of sulphur can be described by coronal 
ionisation equilibrium (Arnaud & Rothenflug 1985^ and that 
the line emission is not in Local thermal equilibrium. The emis- 
sion is plotted on a linear scale (similar to that in figure It 
has been projected onto the sky and convolved with a Gaussian 
ofFWHM = 0'.'l. 

It can be seen that, broadly speaking, there is reasonable 
agreement between the observations and the simulation, at least 
in terms of the morphology close to the head of the bowshock. 
The ring-like features arise from the assumption of cylindrical 
symmetry in the simulation. The emission would, in fact, be 
expected to be a little more 'clumpy' rather than ring-like. 

Finally, the total [SII] luminosity, calculated from the simu- 
lation is 1.1 X 10^^ ergs~^. This is in good agreement with that 
observed (i.e. 4.3 (±3.6) x 10^* erg s^^, see section l3Jl . We 
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Simulation Age 


Corresponding 


Distance 


Bubble 


Longitudinal 


Transverse 


/Age Interval 


Observation 


from source 


Width 


Speed 


Speed 


(years) 


Date/Interval 


(AU) 


(AU) 


(kms ^) 


(kms ^) 


15 


1995 Jan 5 


590 


288 






15-18 


1995 Jan 5 to 1998 Mar 25 


695 


360 


164 


67 


18-19 


1998 Mar 25 to 1999 Feb 3 


722 


386 


148 


63 


19 - 20 


1999 Feb 3 to 2000 Feb 6 


749 


402 


129 


51 


20-21 


2000 Feb 6 to 2001 Feb 10 


775 


420 


115 


44 



Table 4. Projected sizes and speeds of the simulated XZ Tauri outer bowshock as it evolves. Distances are quoted for the latter 
year of the age interval. There is a minimum error of ±6 km s^^ in the simulation speeds. 
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